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Abstract. 

We study the production of heavy quarks, charm at BNL-RHIC (y/s = 200 
GeV/nucleon) and CERN-LHC (y/s = 5.5 TeV/nucleon) and bottom at CERN-LHC 
from heavy ions colliding at relativistic energies. We consider initial fusion of gluons 
(and quark- anti-quark annihilation) , pre-thermal parton interactions and interactions 
in thermalized quark gluon plasma. We also consider free-streaming partons as another 
extreme and compare the results with those from a thermalized plasma of partons. 
The pre-thermal contribution is calculated by considering interaction among partons 
having large transverse momenta (jet-partons) after the initial interaction, and from 
passage of these partons through a thermalized quark gluon plasma. Charm production 
from pre-thermal processes is found to be comparable to that from prompt (initial) 
interactions at LHC. It is suggested that this may have important implications for the 
study of nuclear modification factor, Raa as well as for back-to-back correlation of 
heavy quarks and production of dilcptons having a large mass. 
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Keywords: QGP, heavy quarks, charm, bottom, pre-thermal, back-to-back correlation, 
dileptons. 
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1. Introduction 

Investigation of the properties of quark gluon plasma, a deconfmed strongly interacting 
matter, remains a major activity of present day high energy nuclear physics. This holds 
out a promise of a deeper understanding of the laws of Quantum Chromo-Dynamics 
(QCD) and of early universe which was in the form of quark-gluon plasma. Relativistic 
heavy ion collisions studied at the Brookhaven National Laboratory have produced 
extremely valuable results which confirm elliptic flow of hadrons pQ, suppression of 
hadrons having large transverse momenta or jet-quenching [2], recombination of partons 
as a mechanism of hadronization [3] , radiation of thermal photons jl] , and suppression [5] 
and regeneration [6] of J/ip (see, Ref. [7]). These findings clearly confirm formation 
of quark gluon plasma, which is strongly interacting and has almost a vanishing 
viscosity [H]. 

We are now ready to explore the next and more involved questions about the quark- 
gluon plasma in greater detail and precision. These questions, in no particular order, 
could be the following. When and how quickly does the plasma thermalize? Does 
the time of thermalization depends on the flavors of the partons? What is the flavor 
dependence of energy loss of partons? How does the plasma hadronize? How does 
it evolve in space and time? How does the restoration of the chiral symmetry upon 
formation of the plasma and its breaking upon hadronization affect the evolution of the 
system? How does it break-up? What is the order of phase transition? Is there a tri- 
critical point? Where does it lie? Shall we witness some of the more exotic signatures 
of the formation of the plasma like the local CP violation and chiral condensates? 

A related set of questions concern the production and propagation of the heavy 
quarks in such collisions. Heavy quarks offer some very distinct advantages. They 
are mainly produced from prompt fusion of gluons (gg — > QQ) and quark-antiquark 
annihilation (qq — >• QQ). These processes can be accurately described up to next to 
leading order using perturbative QCD. Their large mass and the necessity to produce 
them as a QQ pair, ensures that their production at later times could be limited. Again, 
due to their very large mass they move slowly, some-what like Gulliver in the Land of 
Lilliput, through the partonic wind of light quarks and gluons. They will lose energy 
as they are buffeted by the partons; through collisions and radiation of gluons. Several 
studies (see, e.g. [9] and references therein) have made detailed evaluations for this 
energy loss. 

Consider b-quarks moving through the plasma. Due to their large mass it is 
expected that while they will lose energy and momentum, the direction of their 
motion may not change substantially due to these interactions which may involve small 
momentum transfers, individually. Their momentum does not undergo a large change as 
they fragment or coalesce with a light quark to form B-mesons. The B-mesons may also 
not change their direction drastically due to their interaction with pions. If true, this 
should make them a valuable probe for the reaction plain dependence of the properties 
of the plasma, especially for non-central collisions. Similar considerations may apply to 
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charm quarks as well, to a great extent. It has also been suggested that the partonic 
flow may sensitively affect the back to back correlations of heavy-quarks [TOl [TT] . 

The modification of the back-to-back correlations reported earlier considers 
production at lowest order [10] and also at NLO [11] , the later contributing substantially 
to the production of charm quarks at LHC. It is quite clear that the heavy-quarks 
produced from splitting of gluons, for example would not move back-to-back. 

If other mechanisms, e.g., a thermal production of heavy-quarks, or a pre- 
equilibrium production of heavy-quarks due to interaction between partons having large 
transverse momenta, or a production due to passage of a parton having a large transverse 
momentum (jet) through the thermalized QGP makes a substantial contribution, then 
these back-to-back correlations will be affected strongly. Of course the so called nuclear 
modification factor Raa will have to be accounted for. It is also expected that the 
relative importance of these contributions will depend on the transverse momentum 
distribution of the heavy quarks, which will introduce further richness in these studies. 
The correlated decay of charm and bottom mesons is also known to lead to a substantial 
contribution to dileptons [12] . which have long been considered a reliable signature of 
quark gluon plasma [T3] . 

In order to put these possibilities to a rigorous test, as a first step we consider 
production of charm quarks at RHIC and LHC energies and of bottom quarks at 
LHC energies due to prompt interactions, thermal productions, and pre-equilibrium 
productions due to interaction of two partonic jets and due to passage of a partonic 
jet through the quark gluon plasma. Our results suggests that there is a substantial 
production of charm quarks at LHC following the initial (prompt) interaction. 

The paper is organized as follows. In the next section we give the formulations of 
prompt, jet-jet, jet-thermal, and thermal interactions leading to production of heavy 
quarks. We shall also consider a scenario, where the initially produced partons undergo 
free-streaming, as an alternative to fully thermalize expansion. The results are discussed 
in Sect. 4. Finally in Sec. 5, we give our conclusions. 



2. Initial Fusion (Prompt Interaction) 



Heavy ion collisions at RHIC (gold on gold) or at LHC (lead on lead) lead to heavy 
quark productions primarily through gluon fusion (gg — > QQ) and also from light quark 
annihilation (qq — > QQ). The flavor excitations for intrinsic heavy quarks (qQ — > qQ 
or gQ — > gQ) is suppressed when next-to leading order processes are considered (see, 
Ref. [H]). 

The differential cross-section for gg — > QQ and qq — > QQ can be written as: 

da \M\ 2 
di 167TS 2 

where the invariant amplitude \M\ 2 is given by 

"(M 2 - if + (M 2 - u) 2 + 2M 2 s 
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Figure 1. The p-r distribution for charm production from initial fusion (solid curve), 
jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and free- 
streaming (solid circles) processes with initial time 0.15 fm/c, in central collision of 
gold nuclei at RHIC at y/s = 200 AGeV. 
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The running coupling constant a s for lowest order is given by 

12tt 
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(33 - 2N f ) ln(Q 2 /A 2 ) ' 

where Nf=3 is the number of flavors and A is the QCD scale parameter used in the 
running coupling constant. 
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Figure 2. The pt distribution for charm production from initial fusion (solid curve), 
jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and free- 
streaming (solid circles) processes with initial time 0.5 fm/c, in central collision of 
gold nuclei at RHIC at y/s — 200 AGeV. The jet-jet contribution in this is thus truly 
pre-thermal. 



The cross-section for the production of heavy quarks from proton-proton collisions 
at leading order [16] is given by 



da 



dyidy 2 d 2 p T 



/i 1 ^ 1 ,Q 2 )/ffe,Q 2 )-* ,l(i ' 



at 



dt 

/a+<y 
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where i and j are the interacting partons and fi and fj are the partonic structure 
functions, and X\ and x 2 are the momentum fractions of the parent nucleons carried 
by the interacting partons. The effect of nuclear shadowing which becomes more 
pronounced at small x has been included using the EKS98 [17] parametrization. We use 
CTEQ5L structure function set for nucleons [18J. The intrinsic transverse momentum 
of the interacting partons is neglected. 

The light quark (u, d, s) and gluon production is similarly calculated with da/di 
taken from Ref. [19] with vanishing quark masses. 

For heavy ion collisions the pt spectrum for heavy quark production is given by 

dN „ da 



d 2 pTdy 



T 



AA 



d 2 pxdy 



(6) 



Heavy Quark Production 



6 




Figure 3. The pr distribution for charm production from initial fusion (solid curve), 
jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and free- 
streaming (solid circles) processes with initial time 0.07 fm/c, in central collision of 
lead nuclei at LHC at ^/s = 5500 AGeV. 



where for central collisions, Tx4=286 fin 2 for Au+Au at RHIC and Taa=292 fm 2 for 
Pb+Pb at LHC. We account for higher order corrections by taking a constant K-factor 
« 2.5 (see, Ref. [9]). 

3. Interaction among Partons 

The initial hard scattering between the partons of the two nuclei will lead to gluons and 
light quarks having large transverse momenta. These will ultimately fragment and lead 
to a stream of hadrons into a narrow cone or jets. We shall continue to call gluons and 
light quarks having large transverse momenta as jet particles. Being copious in number, 
they would interact again and may even approach thermalization. The first collisions 
among two gluonic jets or a quark and anti-quark jet are likely to have sufficient energy 
to produce a pair of heavy quarks. Thus it is felt that jet-jet interaction in relativistic 
heavy ion collision can be an interesting source of production of heavy quarks [20J . 

Several authors have predicted that very large initial temperatures could be attained 
at RHIC [I] and LHC energies [21]. This suggests that even a thermal production of 
heavy quarks through the interaction of thermalized quarks and gluons can take place. 

An interesting new source of high momentum photons [22] and dileptons [23] has 
recently been proposed, which arises from the passage of high momentum quark or gluon 
jets through the thermalized quark gluon plasma. It is of interest to see it whether the 



Heavy Quark Production 



7 



Initial Fusion 




,s 1/2 ={5.5AT^->:; 



V Pb+fj , b(g)LHC : 







1 



2 



3 4 5 



6 7 8 



p T (Gev) 



Figure 4. The pt distribution for charm production from initial fusion (solid curve), 
jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and free- 
streaming (solid circles) processes with initial time 0.5 fm/c, in central collision of 
lead nuclei at LHC at yfs = 5500 AGeV. The jet-jet contribution in this is thus truly 
pre-thermal. 



same process could lead to production of heavy quarks. We shall call this process as jet 
conversion or jet-thermal interaction [24]. 



The general expression for the production of a heavy quark at central rapidity is 
given bv [201 127]: 



where J^P 11 — Pi + P2 — P ~ p'-> Pi an d P2 are the four-momenta of the incoming partons 
and p and p' are the same for outgoing heavy quark and anti-quark. F(x,p, t) gives the 
phase space distribution function for the incoming partons. 

Before deriving results for different processes by choosing appropriate phase space 
distributions etc., we can perform the following simplifications: 

Writing d 3 pi/iJi = p^dpT^idyi and integrating over d 3 p', we get for y=0 (p z =0, 




x* 4 (Ep m ) /e' 

x |M| 2 F(f,pl,t)F(x 2 ,p 2 ,t) 



(7) 
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Figure 5. The pt distribution for bottom production from initial fusion (solid 
curve), jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and 
free-streaming (solid circles) processes with initial time 0.07 fm/c, in central collision 
of lead nuclei at LHC at V« = 5500 AGeV. 



and Pt=p)' 

I (2>") M F(i,n,t)F(l,n,t) 

= J dyi dy 2 &<\>\ dfa PT2dp T 2 Pridpri ^ \M\ 2 
xF(x,pri,<l>i,yi,t)F(x, p T 2, 02, 1/2, *) , 

(8) 

where 

6 (£E) S(pti-Pti,o) 



E' \p T2 (cosh(y 1 - y 2 ) - cos(0i - fe)) - (E cosh y x -p cos fa)] 

and we have 



(9) 



PT2 (E cosh y 2 -p cos 2 ) 

T1 '° [pT2(cosh(yi - j/ 2 ) - cos(0i - 2 )) - (E cosh j/ x - p COS (pi)] 
Now we proceed to evaluate individual contributions. 



3.1. Jet- Jet interaction of partons 

In order to estimate the contribution of the jet-jet interaction to the production of heavy 
quarks, we approximate [20] the phase-space distribution of the gluon, quark, or anti- 
quark jets produced in initial (prompt) scattering of the partons in a central collision 
as: 

F(x,p,t) = f jet (x,p,t) , (11) 
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Figure 6. The pt distribution for bottom production from initial fusion (solid curve), 
jet-jet (dashed curve), jet-thermal (dash-dotted curve), thermal (stars), and free- 
streaming (solid circles) processes with initial time 0.5 fm/c, in central collision of 
lead nuclei at LHC at s/s = 5500 AGeV. The jet-jet contribution in this is thus truly 
pre-thermal. 



fi*& a *) = .loL 6 (y - v) Q ( r / - T ) ( r - ^ > ( 12 ) 



where [20] 

(2tt) 3 dNi 
giTirfirPT dyd 2 p T 

with pr > 2 GeV. Here 5(y — if) denotes the Bjorken correlation for space-time and 
energy- momentum rapidities and V stands for quarks, anti-quarks, or gluons. The 
degeneracy of quarks and gluons is given by g g / q such that g g = 8 x 2 and g q = 3 x 2 |25j. 
Further is the transverse radius of the nucleus and dNi/ 'd^p^dy is the transverse 
momentum distribution of partons for pr > 2 GeV. We neglect the dependence of this 
distribution on the momentum rapidity [20] as we are calculating the results for heavy 
quarks at y=0, when only very small values of yi and j2 contribute, and the rapidity 
dependence is marginal. For results at y ^ 0, appropriate distributions will need to be 
used. 

The momentum space distribution of the jets at RHIC and LHC are taken from 
parametrization [23] given earlier, where the jet distributions were calculated in LO- 
pQCD with a K-factor (« 2.5), CTEQ5L structure functions and EKS98 shadowing 
functions, which we have used here for calculation of initial production of heavy quarks. 
Thus we have, 



dN da iet 
Taa- 



dyd 2 p T d?p T dy 



r _r K ^hw (13) 
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and 



1 dN 



Qi d?p T dy 



(14) 



y=0 



where K, C, B, and j3 are taken from reference 

We neglect transverse expansion of the system, which should be valid at early times 
when most of the heavy quarks are produced. Now taking 

d 4 x = t drr dr dr) d<p r , (15) 

we can perform the integration over r, <fr r and r. 

Thus the pr distribution of open heavy quark production from jet-jet interaction 
can be written as 

J 3 N Hv/Ti) 



E 



d 3 p 



16(27T) 2 (7Ti?2) 



drj dp< 



T2 



X 



Pt2(1 — cos(0i — 02 )) — (Ecoshr) — pcos0i} 



1 



y9g M et (Pn,0) ^2 )\ M \gg^QQ 

+g 2 q N f hl t (p Tlfi )h] ct (p T2 ) \M\\^ Q 
3 is the number of quark flavors [26] . 



(16) 



where Nf 

While writing the above, we have used the Bjorken correlations 5(y\ — rj) and 
S(y 2 - rj) (see EqJTSJ. With this, the Eqs.® and (JTOD reduce to, 
S(£E) 5(p T1 -p Tlfi ) 



and 



E' 



Pt2(1 — cos(0i — 2 )) — (-Ecoshr^ — pcos0i) 
PT2 (-E 1 cosh rj — p cos fa) 



;i7) 



PT2(1 — cos(0i — 4>2)) — (E cosh T) — pcos0i) 
and thus Eq.® reduces to Eq. ffTBT) . above. Numerical integration of the EqJTBI gives 
pre-thermal heavy quark production. We shall see that the major contribution to heavy 
quark production comes from gluon fusion. 

We add that Levai et al. [27J calculated the so-called pre-equilibrium (jet-jet) 
contribution by assuming that this mechanism operates during r e [0.1 - 0.5] fm/c. 
We shall explore some other options as well. 

We re-iterate that the above expressions do not account for energy loss of the 
energetic gluons and quarks as they traverse the plasma and thus they provide the 
upper limit of the heavy quark productions. 



3.2. Interaction among thermally and chemically equilibrated partons 

If multiple scatterings occur in a rapid succession, QGP may reach thermal equilibrium 
quite early and follow hydrodynamics evolution till it hadronizes. During this phase of 
evolution, heavy quarks may be produced provided the initial temperature is high [28J. 
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We can estimate the initial temperature by assuming Bjorken hydrodynamics 
which relates it to the particle rapidity density by 

2vr 4 dN _ , . 

4aT 3 r (19) 



45C(3)7r J Rf, dy 

where dN/dy is the particle rapidity density, a=42.257r 2 /90 for massless light quarks 
and gluons and Rt = 1.2v4 1//3 . We take particle rapidity density as 1260 estimated 
experimentally [30J at RHIC and assumed [21] that the particle rapidity density at LHC 
is about 5625. Some recent works suggest a smaller value for dN/dy ~ 3000-3500 at 
LHC [31] from considerations of parton saturation. However, larger values have also 
been suggested j32|. The initial experimental results from pp collisions at LHC at 0.9 
TeV, 2.36, and 7 TeV already show an increase in the particle rapidity densities, which 
is steeper than expected [33]. In any case, results for any other rapidity density can be 
easily obtained. 

The time evolution of the temperature of thermalized QGP for a boost-invariant 
longitudinal expansion is governed by [29J: 

T 3 r = T 3 t = const. (20) 

Assuming a rapid thermalization and chemical equilibration, we get the lowest estimate 
of r from the above by taking r ~ 1/3T [21J. This provides that r ~ 0.15 fm/c 
at RHIC. Such a small value for r is supported by the single photon data [I]. As an 
alternative, we also consider a much larger time of thermalization, tq ~ 0.5 fm/c, with 
T calculated from Eq. [191 We add that several studies, especially the ones related to 
the flow of hadrons tend to use a larger value of r ~ 0.6 fm/c p], even though more 
recent studies find only a weak dependence of r on the flow for hadrons [31]. For LHC, 
we have similarly assumed r ~ 0.07 fm/c and 0.5 fm/c. 

Taking the critical temperature is taken to be 0.170 GeV, we see that the end of 
the QGP phase occurs at: 

r f = (const. /0.170 3 ) (21) 

and then the thermal production mechanism would operate during tq to t/. 

We take the phase space distribution for the thermalized quarks and gluons as, 

fthiPr, y> V) = exp [ Pt cosh(?/ - rj)/T] . (22) 
Thus the transverse momentum distribution of thermally produced charm given by 



d?p 



16(2vr) 



J r dr dr] dprp 2 d(j)\ d(j) 2 dy\ dy 2 
(PT2 Pti,o) 



[p T2 (cosh(?/i - y 2 ) - cos(0i - 2 )) - {Ecos\iy l -pcos0i)] 

x [fth(pTi,a,yi,v)fth(pT2iy2,v)) x 



9g 



M 99 ^QQ 



+ 9,N } 



M qq^QQ 



(23) 
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Just as in pre-thermal case, we calculate the final charm production at y=0 or 
central rapidity and thus p z = and px = P- We now have the kinematical constraint; 

S{J2E) 5(pti-Pti,o) 



E> 



and 



PTlfi 



\p T2 (cosh(y 1 - y 2 ) - cos(0i - fa)) - (Ecosh?/! -pcos^i)] 
(PT2{E cosh y 2 -pcos0 2 )) 



(24) 



(25) 

[p T2 (cosh(y 1 - y 2 ) - cos(0i - 2 )) - (Ecosh?/! -pcos^i)] 

Numerical integration of the Eq. for different initial conditions gives us the 
contribution from thermalized QGP. 

3.3. Interactions of free- streaming partons 

As an extreme, we consider free-streaming partons, as a model of evolution of the 
system of deconfined quarks and gluons, which completely relaxes the condition of 
thermalization. The initial distribution at t = tq and z = is obtained by assuming 
maximum entropy [21], so that 



dN 



exp(-— ) 

-L n 



d 3 pd 3 x 

and the condition that needs to be satisfied is 
^ df(x,p) 
dx^ 



pr — : — — = . 
We assume boost invariance along the z-axis with 

f(p,x) = f(p T ,p z t-Ez) . 
The solution which satisfies the differential Eq. |27] is 



f(p, x) = exp 



p\ + (p z t - Ezf/ri 



(26) 

(27) 
(28) 

(29) 



Eq. 



Now using 

Pz 

Z 

becomes 



p t sinh y, E = p t cosh y, 
rsinhr/, t = r cosh n , 



f(PT,V,y) = exp 



Pxyl + T 2 sinh (y — rj)/n 



(30) 



(31) 



Thus the final integration to calculate pr distribution for heavy quark production 
from free streaming partons is given by 
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E 



d 3 N 



d 3 p 



TlRl 



y=0 



16(2 



7T 



J rdr drj dpT2 d<j)i d<\>2 dyi dy 2 

(PT2 PTlfl) 



[pT2{cosh(y 1 - y 2 ) - cos(0i - 2 )) - {Ecoshy 1 -pcos0i)] 



X [f(PTl,0,Vl,Vl) f{PT2,V2,V2)} X 



o9g 



M. 



gg->QQ 



+ 9 2 g N f 



M, 



qq^fQQ 



(32) 



The initial conditions for the free-streaming case are taken to be same as that for the 
thermal production, whereas the final time is taken as Rt/c, the transverse radius of 
the nuclei, after which the system would surely expand rapidly along the transverse 
direction as well and disintegrate. 



3.4- Charm production from the passage of jets through thermal medium 

Now we discuss the production of heavy quarks by passage of light quark and gluonic 
jets through thermalized QGP. In order to proceed, we can differentiate the phase space 
distribution into jet partons and thermalized partons. The phase space distribution for 
equilibrated medium is given by 

fth = exp [-p T cosh(?/ -T))/T] . (33) 

The jet distribution is given by Eq. [11] and [12] for px > 2 GeV. We have already 
discussed Tj or r , which gives the start of the time from when we consider the system to 
be in the form of QGP. We define r d as the time which a jet takes to reach the surface 
of the quark gluon plasma. Consider a jet formed at r*with velocity v which travels to 
the surface of plasma. The distance, d, covered in this process is given by, 

d = —r cos + \J R\ — r 2 sin 2 , (34) 

where = cos -1 (t> • f), and Rt is the radius of the system. A massless quark or a gluon 
would take a time 

r d = d/c (35) 

for this journey. Considering that QGP would cool down to the critical temperature 
by r = Tf, the time spent by the jet in the plasma would be r max = min [rf,T d ] (see 
Ref. [23]), and the jet-thermal production mechanism would be in operation during r e 

i > ''"max] • 

The thermalization time is taken as either 0.15 fm/c or 0.5 fm/c at RHIC and as 
either 0.07 fm/c or 0.5 fm/c at LHC as before. 

The final result for jet-thermal interaction is given by 
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d 3 N 



d 3 p 



y=0 



16(27T) 4 7T^ 



J dr rdr drj d(j)\ d(f>2 dy\ dpx2 



(Pr2(cosh(?/ 1 — 7}) — cos(0i — 02 )) — E cosh yi + p z sinh y x + p T cos < 



xfth{PTi,o,Vi,v) 



9 2 q N f hU PT2 ) \M\'n + -g 2 g hU PT2 ) \M\ 



1 99- 



(36) 



which we evaluate numerically. 

We add that these results do not account for the energy loss suffered by the high 
energy quark of gluon as it traverses the plasma, and during which it may even change 
its flavour. This would necessitate a treatment along the lines of Ref. [35] for production 
of photons from a similar process. 



4. Results 

Now we discuss our results. As a first step we show the results for charm production at 
RHIC and LHC. 



4.1. Charm at RHIC 

In Fig. [H we plot the results at RHIC for charm production from prompt interactions 
(initial fusion), jet-jet interactions of partons, thermal production and from the passage 
of high momentum jets through thermalized quark-gluon plasma. We also give the free- 
streaming results to compare with thermal production. We consider central collision of 
gold nuclei at V« = 200 AGeV. 

We see that the contribution of the initial fusion dominates at all pt- Considering 
initial formation time tq as 0.15 fm/c, we can consider two extremes for the jet-jet 
interactions as indicated earlier. For one extreme we consider that jet-jet interaction 
continues till jets reach the surface of the system (77 ~ Rt/c). The other extreme could 
be to assume that it operates only till the time of thermalization, if it is large (t/ ~ 0.5 
fm/c), so that it is considered as the pre-equilibrium contribution (27]. In any case the 
time integration for this contribution reduces to In (tj/tj), (see Eq. [16]) and thus one can 
easily obtain this results for any choice of initial conditions. 

Our results for this contribution at RHIC are similar in magnitude and form to 
those reported by Lin and Gyulassy [20J for RHIC energies. Recall that we consider 
only partons having px > 2 GeV as constituting the jets. The jet-thermal contribution 
is seen to be comparable to the contribution of the jet-jet interaction. 

We find that the contribution of thermal production at large pt is rather small. 
However it is larger by a factor of about 3 at lower p?, compared to the jet-jet and 
jet-thermal contributions. This has its origin in large initial temperature for low r 
value assumed here and exclusion of partons having px < 2 GeV in the jet distribution. 
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The contribution obtained by assuming free-streaming partons having same initial r 
and operating till tj = Rt/c is quite similar to thermal contribution. We also note 
that both thermal and free-streaming contributions drop a lot more rapidly for larger 
Pt compared to contributions initiated by prompt and jet interactions. The latter 
contribute at a level of about 4-5 % for large pt as compared to the initial fusion for 
tq = 0.15 fm/c and T = 447 MeV. We add that a much larger initial temperature of 
700 MeV and K=4 has been assumed in the calculations by Liu and Fries [24J for the 
jet conversion in thermalized medium. 

In Figj2j we give results obtained by using the formation time of the plasma as 0.5 
fm/c, and further assuming that the jet-jet interaction is in operation during r « 0.15 
fm/c to 0.5 fm/c, so that it can be considered as a pre-equilibrium contribution (27]. In 
any case, we find that initial fusion gives the largest contribution to charm production, 
followed by jet-jet and jet-thermal interactions. At high pr, jet-jet as well as jet-thermal 
interactions give a contribution of about 1% of the initial fusion. Now the thermal and 
free-streaming contribution though remaining similar in magnitude, are much smaller 
and also fall more rapidly with p?. They give a contribution at a level of 1% of the 
initial fusion contribution at pr < 3 GeV. 

We realize that even though small in magnitude, the jet-jet and jet-thermal 
contributions may still give a discernible feature to the back-to-back correlation of 
charm-quarks at RHIC. The correlation of charm quarks from these should be distinct 
from initial fusion (back-to-back for LO contribution and forward peaked for NLO 
contribution). This aspect is under study. 

We have not given our results for bottom production at RHIC energies here as we 
have found that bottom production from the jet-jet mechanism is less than 3 orders 
of magnitude, the thermal mechanism is less than 6 orders of magnitude and the jet- 
thermal mechanism is less than 4 orders of magnitude of that obtained from initial 
(prompt) interaction [9]. We do add, however, that the decay of bottom into charm 
quarks adds an interesting richness to the study of charm production as well as charm- 
correlation. 

4-2. Charm and bottom at LHC 

Next we discuss our results for charm and bottom production at LHC. 

We consider central collision of lead nuclei at = 5500 AGeV. Other initial 
conditions have already been discussed. We find that (see Figj3]) the charm production 
from initial fusion is about a factor of 10 or more than that at RHIC, and of-course its 
fall with pt is considerably slower, as one would expect. 

We plot our results at LHC, taking time of formation of QGP to be 0.07 fm/c. We 
have assumed the initial time of jet-jet interaction to be from < 0.1 fm/c until the 
jets reach the surface of the system, Tf rs Rt/c as one of the extremes. We find jet-jet 
contribution even exceeds the initial fusion for px > 2 GeV. Of-course since we include 
jets having p? > 2 GeV, this contribution decreases for lower p? ■ We shall return to 



Heavy Quark Production 



16 



this. The thermal production of the charm, taking advantage of the initial temperature 
when Tj ~ 0.07 fm/c, is almost 40% of the initial fusion contribution at low pt, and 
drops to about 5-10% of the same at pt ~ 5 GeV. The jet-thermal contribution is also 
seen to be considerable at large p?, though smaller than jet-jet contribution. It is seen 
to be comparable to the thermal contribution at lower p?. Of-course this has its origin 
in the large initial temperature. 

Next in FigJH we plot our results at LHC for another set of initial conditions. 
Taking time of equilibration to be 0.5 fm/c, we have pre-thermal charm production 
resulting from jet-jet interactions starting at Tj = 0.07 fm/c and operating till the time 
of thermalization. Jet-jet interaction gives a contribution which is comparable to that 
from initial fusion. Once again the decrease at lower p? occurs due to the exclusion 
of jets having px < 2 GeV. We also note that for px ~ 8 GeV, the jet-jet interaction 
contributes at a level of 10% of the initial fusion. The jet-thermal contribution on the 
other hand starts at about 20 % of the initial fusion contribution at the lowest pt and 
rises to about one-third of initial contribution for larger p?- The thermal and the free- 
streaming contributions are seen to be quite small for larger p?, though around p? ~ 
GeV, they contribute at a level of 10% of the initial fusion. 

The thermal and the free-streaming contributions shown in Figs. |3] and H] deserve 
more attention. We see that the two contributions for larger initial temperatures at 
LHC differ by a factor of about 2 or more, while for the smaller initial temperature 
at RHIC they are of similar magnitude. This we feel has its origin in the large initial 
temperature which enhances the phase-space contribution to the thermal production of 
charm quarks. In fact in Figs. [5] and [6j where we plot the contributions for the bottom 
quarks production at LHC, this difference again shows up. 

Coming back to our earlier observation about the contribution of these processes 
to the back-to-back correlation of charm quarks, we realize that the large contribution 
from processes like jet-jet, jet-thermal and even thermal production of charm quarks at 
LHC is likely to drastically alter the conclusions about this which were arrived at from 
studies [TOT ITT] which invoked only the initial fusion processes. 

As indicated earlier, we give our results (see Figs. |5] and [6]) for the production of 
bottom quarks at LHC for two initial conditions discussed earlier. Of-course we find 
that the thermal production of bottom quarks is quite negligible at LHC. The jet-jet 
contribution is seen to be at a level of 2-6% of the initial fusion. The jet-thermal 
contribution is also negligible. 

Thus we feel that a high statistics data may still be able to discern the contribution 
of processes other than initial fusion to the back-to-back correlation. We shall report 
this in our future publications. 

5. Summary and Conclusions 

We have reported results on charm production at RHIC and LHC and bottom 
production at LHC from initial fusion and multiple scattering processes treated as jet- 
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jet interaction and thermal production and passage of jets through QGP. Two different 
initial conditions, one with early thermalization and other with thermalization at n ~ 
0.5 fm/c have been used. Substantial production of charm, specially at LHC, is seen in 
addition to the production due to initial fusion. 

The jet-jet contribution to charm quarks at LHC even exceeds the initial fusion 
contribution for intermediate for suitable initial conditions. We have argued that 
since back-to-back correlation for charm quarks from the processes under consideration 
could be different from those from initial fusion, the results for that, specially at 
LHC would be considerably affected. This may even be discernible for back-to-back 
correlation of bottom quarks at LHC. As indicated earlier, these results will be published 
shortly. 

One could think of several improvements. The obvious one would be to include 
energy loss suffered by the interacting gluons and quarks before fusion and those by 
the heavy quarks after production. We feel that the inclusion of energy loss before 
fusion may have smaller effect on over-all production of the heavy quarks as these 
would be limited to the earliest times when the momenta and the temperature are 
still very large. The energy loss suffered by the heavy quarks after the production can 
not be neglected (see e.g., Ref. [9]). This will alter the px distribution of the heavy 
quarks. However, as it will affect the heavy quarks produced by all the processes, 
their relative contributions will remain largely unaffected. Relaxing the condition of 
chemically equilibrated plasma [27] will also be useful. It should be of interest to, 
alternatively, estimate the prompt charm production using the colour glass condensate 
model which predicts different scalings for different rapidities [36] . 

A more complete calculation using parton cascade model [37] will be reported 
shortly (see also, Ref. [38]). 

We conclude that production of charm quarks at LHC and even at RHIC from 
processes other than initial fusion can be large and can play a significant role in our 
study of back-to-back correlation. This will have important implications for the study 
of the nuclear modification factor Raa as well as large mass dileptons having their origin 
in the correlated charm decay [12, 28J. 
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